LARGE CLIQUES IN A POWER-LAW RANDOM GRAPH 



SVANTE JANSON, TOMASZ LUCZAK, AND ILKKA NORROS 

Abstract. We study the size of the largest clique o»(G(n, a)) in a ran- 
dom graph G(n, a) on n vertices which has power-law degree distribution 
with exponent a. We show that for 'flat' degree sequences with a > 2 
whp the largest clique in G(n, a) is of a constant size, while for the 
heavy tail distribution, when < a < 2, uj(G(n, a)) grows as a power 
of n. Moreover, we show that a natural simple algorithm whp finds in 
G(n, a) a large clique of size (1 + o(l))ui(G(n, a)) in polynomial time. 



1. Introduction 

Random graphs with finite density and power-law degree sequence have 
attracted much attention for the last few years (e.g. see Durrett [§[ and 
the references therein). Several models for such graphs has been proposed; 
in this paper we concentrate on a Poissonian model G(n, a) in which the 
number of vertices of degree at least i decreases roughly as ni~ a (for a precise 
definition of the model see Section [2] below). 

We show (Theorem [T]) that there is a major difference in the size of the 
largest clique u(G(n, a)) between the cases a < 2 and a > 2 with an 
intermediate result for a = 2. In the 'light tail case', when a > 2 (this is 
when the asymptotic degree distribution has a finite second moment), the 
size of the largest clique is either two or three, i.e., it is almost the same 
as in the standard binomial model of random graph G(n,p) in which the 
expected average degree is a constant. As opposite to that, in the 'heavy 
tail case', when < a < 2, ui(G(n,a)) grows roughly as n 1- "/ 2 . In the 
critical case when a = 2 we have to(G(n, a)) = O p (l), but the probability 
that G(n, a) > k is bounded away from zero for every k. We also show 
(Corollary [3]) that in each of the above cases there exists a simple algorithm 
which whp finds in G(n, a) a clique of size (l — o(l))uj(G(n,a)). This is quite 
different from the binomial case, where it is widely believed that finding large 
clique is hard (see for instance Frieze and McDiarmid 

Similar but less precise results have been obtained by Bainconi and Mar- 
sili [H; 0] for a slightly different model (see Section 16.61 below) . 
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2. The model and the results 



The model we study is a version of the conditionally Poissonian random 
graph studied by Norros and Reittu 15] (see also Chung and Lu [7] for a 



related model). For a > 1 it is also an example of the 'rank 1 case' of 
inhomogeneous random graph studied by Bollobas, Janson, and Riordan 0, 
Section 16.4]. 

In order to define our model consider a set of n vertices (for convenience 
labelled 1, .. . ,n). We first assign a capacity or weight Wi to each vertex i. 
For definiteness and simplicity, we assume that these are i.i.d. random vari- 
ables with a distribution with a power-law tail 

F(W > x) = ax~ a , x > x , (2.1) 

for some constants a > and a > 0, and some xo > (here and below 
we write W for any of the Wi when the index does not matter). Thus, 
for example, W could have a Pareto distribution, when xq = a x l a and 
F(W > x) = 1 for x < xq, but the distribution could be arbitrarily modified 
for small x. We denote the largest weight by W mas = max^Wj. Observe 
that (12.11) implies 

HW max > tn 1/a ) < n¥(W > tn 1/a ) = 0(r a ). (2.2) 

Note also that E W 13 < oo if and only if a > (3; in particular, for the 'heavy 
tail case' when a < 2 we have EVF 2 = oo. 

Now, conditionally given the weights {V^}™, we join each pair of 
vertices by Ey parallel edges, where the numbers are independent Pois- 
son distributed random numbers with means 

EE ij = X ij = b ] ^-, (2.3) 
n 

where b > is another constant. We denote the resulting random (multi) graph 
by G(n, a). For our purposes, parallel edges can be merged into a single edge, 
so we may alternatively define G(n, a) as the random simple graph where 
vertices i and j are joined by an edge with probability 

Pij = 1 - exp(-Ay), (2.4) 

independently for all pairs with 1 < i < j < n. 

Then our main result can be stated as follows. Let us recall that an event 
holds with high probability (whp), if it holds with probability tending to 1 
as n — > oo. We also use o p and O p in the standard sense (see, for instance, 
Janson, Luczak, Rucihski [HI]). 

Theorem 1. (i) If < a < 2, then 

u(G(n,a)) = (c + 0p (l))n 1 - Q / 2 (logn)- a / 2 , 

where 

c = ab a/2 {l-a/2)- a ' 2 . (2.5) 
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(ii) If a = 2, then uj{G{n,a)) = O p (l); that is, for every e > there 
exists a constant C £ such that P(uj(G(n, a)) > C £ ) < e for every n. However, 
there is no fixed finite bound C such that u>(G(n,a)) < C whp. 

(Hi) If a > 2, then Lo(G(n,a)) £ {2,3} whp. Moreover, the probabilities 
of each of the events u>(G(n,a)) = 2 and u>(G(n,a)) = 3 tend to positive 
limits, given by (|5.10p . 

A question which naturally emerges when studying the size of the largest 
clique in a class of graphs is whether one can find a large clique in such 
graph in a polynomial time. By Theorem whp one can find co(G(n,a)) 
in a polynomial time for a > 2, and, with some extra effort, the same can 
be accomplished for a = 2 (see Corollary [3]) . Thus let us concentrate for 
the case a < 2, when the large clique is of polynomial size. Let us suppose 
that we know the vertex weights Wi defined in Section [2] and, for simplicity, 
that these are distinct (otherwise we resolve ties randomly; we omit the 
details). Since vertices with larger weights tend to have higher degrees, 
they are more likely to be in a large clique, so it is natural to try to find 
a large clique by looking at the vertices with largest weights. One simple 
way is the greedy algorithm which checks the vertices in order of decreasing 
weights and selects every vertex that is joined to every other vertex already 
selected. This evidently yields a clique, which we call the greedy clique and 
denote by /C gr . Thus 

/C gr = {i : i ~ j for all j with Wj > Wi and j E /C gr }. 

A simplified algorithm is to select every vertex that is joined to every 
vertex with higher weight, regardless of whether these already are selected 



or not. This gives the quasi top clique studied by Norros 14i |. which we 
denote by /C qt . Thus 

^Cqt = {i ■ i ~ j for all j with Wj > Wi}. 

Obviously, /C qt C /C gr . The difference between the two cliques is that if we, 
while checking vertices in order of decreasing weights, reject a vertex, then 
that vertex is ignored for future tests when constructing /C gr , but not for 
/Cqt- A more drastic approach is to stop at the first failure; we define the 
full top clique /Cf t as the result, i.e. 

/Cft = {i : j ~ k for all distinct j, k with Wj,Wk > Wi}. 

Thus /Cft is the largest clique consisting of all vertices with weights in some 
interval [w, oo). Clearly, /Cf t C jC qt C fC gT . Finally, by /C max we denote the 
largest clique (chosen at random, say, if there is a tie). Thus 

|/Cft| < |/C q t| < |/C gT | < |/C max | = tj(G(n, a)). (2.6) 



The following theorem shows that the last two inequalities in (|2.6p are as- 
ymptotic equalities, but not the first one. Here we use — — > for convergence 
in probability, and all unspecified limits are as n — ► oo. 
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Theorem 2. If < a < 2, then G(n,a), /C gr and /C qt both have size 
(1 + o p (l))uj(G(n, a)); in other words 

\Kgr\/\K max \ 1 and |/C qt |/|/C ma x| 1- 
On the other hand, 

|ACf t |/|/C max | -^-> 2~ a l 2 . 

Thus, whp /Cg r and /C q t almost attain the maximum size of a clique, while 
JCft falls short by a constant factor. As a simple corollary of the above result 
one can get the following. 

Corollary 3. For every a > there exists an algorithm which whp finds in 
G(n,a) a clique of size (1 + o(l))a;(G(n, a)) in a polynomial time. 

3. The proof for the case a < 2 (no second moment) 

We begin with a simple lemma giving an upper bound for the clique 
number of the Erdos-Renyi random graph G(n, y) (for much more precise 
results see, for instance, Janson, Luczak, Rucihski [13l]). 

Lemma 4. For any p = p{n), whp 

mi \\ s 2 log re 

u(G(n,p)) < . 

1 — p 

Proof. Denote by the number of cliques of order k in G(n,p). For the 
expected number of such cliques we have 

'k\ / n.p , n /„\ k 

~ v k 

If we set k > [2 log(n) /(l — p)\ , then 

p (*-i)/2 = (!_ {1 _ p) f-m < e -(i- p) ( fe -i)/2 < e/n> 

Consequently, we arrive at 

F(u(G(n,p)) >k)= F(X k > 1) < EX k < {^f -> 0, 
since > |_21ognJ — ► 00. □ 

Proof of Theorems a^rf [H For s > 0, let us partition the vertex set 
V = {1, . . . , n} of G n = G(n, a) into 



iE^ = c;)^<(yp (fc - i)/2 



V, = {i : Wi < s\/n log n} and y, + = {i : Wj > s-\/n logn}; 

we may think of elements of V~ and V s + as 'light' and 'heavy' vertices, 
respectively. By (|2.ip . 

E|F+| = nP(W > s^nlogn) = as~ a n 1 ^ 2 log- a/2 n. (3.1) 

Moreover, |V^ + | ~ Bin(n,P(VF > s\/n logn)), and Chebyshev's inequality 
(or the sharper Chernoff bounds (l3l . Section 2.1]) implies that whp 

\V+\ = (l + o(l)) E\V s + \ = {1 + o{l))as- a n l - a/2 \og- a/2 n. (3.2) 
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We now condition on the sequence of weights {Wjt}. We will repeatedly 
use the fact that if i and j are vertices with weights W% = x\Jn logn and 
Wj = y\Jn logn, then by (|2.3p - (|2.4p . A^- = bxy logn and 

Pij = 1 - e ~ x ^ = 1 - e - bxylosn = 1 - n- bxy . (3.3) 

In particular, 

Pij <l-n- bs \ ifi,jeV~, (3.4) 

Pij >l-rr bs \ iii,jeV s + . (3.5) 

Consider, still conditioning on {W&}, for an s that will be chosen later, the 
induced subgraph G n [V~] of G(n, a) with vertex set V~. This graph has at 
most n vertices and, by (|3.4p . all edge probabilities are at most 1 — n~ bs , so 
we may regard G n [V s ~] as a subgraph of G(n,p) with p = 1 — n~ bs . Hence, 
Lemma H] implies that whp 

u{G n [V-])< 2 -^ = 2n bs2 logn. (3.6) 

If K, is any clique in G(n,a), then K, n V~ is a clique in G n [V r s _ ], and 
thus \ICnV~l < u(G n [V s -])) further, trivially, \)CnV s + \ < \V+\. Hence, 
\JC\ < u(G n [V-]) + \V+\, and thus 

u;(G(n,a))<u;(G n [V-]) + \V s + \. (3.7) 

We choose, for a given e > 0, s = (1 — e)6 _1 / 2 (l — a/2) 1 / 2 so that the 
exponents of n in (13.6H and (j3.2|) almost match; we then obtain from (|3.7p . 
(1331) . and QMSl, that whp 

w(G(n,a)) < (l + o{l))as- a n l - a / 2 log- a / 2 n 

= (1 + o(l))(l - e)" ^ 1 - / 2 log- Q/2 n, (3.8) 

with c defined as in (|2.5p . 

To obtain a matching lower bound, we consider the quasi top clique /C qt . 
Let, again, s be fixed and condition on the weights If i,j £ Vj~, 

then by (|3.5p . the probability that i is not joined to j is less than n~ bs . 
Hence, conditioned on the weights {W&}, the probability that a given vertex 
i G Vg + is not joined to every other j £ is at most |V s + |n _b ' 52 , which 
by (|3.2p whp is at most 2as~ a n l ~ a l 2 ~ bs log~ a ^ 2 n. We now choose s = 
(1 + e)b~ l l 2 (l — a/2) 1 / 2 with e > 0. Then, for some constant C < do, whp 

P(i £ /C qt | {W k }) < Cn- 2 ^- a W 

and thus 

H\V S + \ /C qt | I < Cn-^-^lV+l. 

Hence, by Markov's inequality, whp 

|^ + \£ q t| < Cn~ £ ^- a ^\V s + \. 
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Thus, using (j3.2j) . whp 

u(G(n,a)) > |/C gr | > |/C qt | > \V+\ - \V+\lC qt \ > (1 - o(l))\V+\ 

> (1 - o{l))as- a n 1 ~ a / 2 log~ a/2 n 
= (1 + o(l))(l + e^cn 1 -"/ 2 log- a/2 n. (3.9) 



Since e > is arbitrary, (|3.8p and (|3.9[) imply Theorem [ J^i)| and the first 
part of Theorem [2j 

In order to complete the proof of Theorem [2j it remains to consider /Cft. 
Define G n as the complement of G(n, a). Then, using (|3.5p and conditioned 
on {Wj}, we infer that the expected number of edges of G n with both end- 
points in is at most n~ bs |V S + | 2 . If we choose s = 6 _1//2 (2 — a) 1 / 2 , then 
(|3.2p implies that this is whp o(l); hence whp K+ contains no edges of G n , 
i.e., /Cft D V+. 

On the other hand, let < e < 1/2 and define, still with s = b~ 1 / 2 (2 — 
a) 1//2 , V = y^_2 £ ) s H V"(7-e)s - Then, conditioned on {Wi} ; the probability 



of having no edges of G n \\\V is, by ([3 

[I ftj < (1 - n-^ 1 - 6 ) 3 - 8 )^') < exp(-n-( 1 - £ ) 2 ( 2 - Q )(|y'| - l) 2 /2). 
i,jeV 

(3.10) 

By (j32J), whp 



\v'\ - 1 = | - iy+ i-i 

1 1 1 (1— 2e)s 1 1 (1— e)s< 

= (1 + o(l))a((l - 2e)~ a - (1 - e)" )^ /! 1 - / 2 log~ a/2 n, 
and it follows from (|3.10p that 

p(/c ft ^ v+_ 2£)s ) < n ^ -> °- 

Hence, whp /C ft C V£|_ 2e y 

We have shown that, for any e € (0, 1/2), whp C /Cf t C Ki_ 2e ) s j an d 
it follows by <g£2j) and ([23]) (by letting e -> 0) that whp 

|/C ft | = (1 + o(l))\V+\ = (1 + o(l))a S - a n 1 -"/ 2 log""/ 2 n 

= (1 + (l))2~ a / 2 cn 1 - a / 2 log""/ 2 n = (1 + o(l))2-"/ 2 ^(G(n, a)), 
where the last equality follows from Theorem [TJ □ 

4. The case a = 2 (still no second moment) 

Proof of Theorem\ ^\\)\ and Corollary^ Given the weights Wj, the proba- 
bility that four vertices i,j,k,l form a clique is, by (|2.4p and (|2.3|) . 

PijPikPilPjkPjlPkl < KjKkKl^jk^jl^kl — b . 
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Thus, if X m is the number of cliques of size m in G(n, a), then the conditional 
expectation of is 



E(^4 | {Wi}^) < 6 6 n~ 6 WfWfWlWf < b 6 (n~ 3 / 2 Wf J . (4.1) 

i<j<k<l i 

To show that the number of such quadruples is bounded in probability, 
we shall calculate a truncated expectation of Y^i^f- Using (|2.ip . for any 
constant A > 0, we get 

e(J2 W ?; ^max < ^™ 1/2 ) < E(j]min(Wi, An 1 ' 2 ) 3 ) 

= nEmin(PF, An 1/2 ) 3 

= n / 3x 2 P(W > x) dx 
Jo 

= 0(nAn 1/2 ), (4.2) 

and thus, using (|2.2p and Markov's inequality, for every i > and some 
constant C independent of A, t and n, we arrive at 



i 

< t- 1 E(V 3 / 2 Y Wf- W max < An 1 ^ + P(iy max > An 1 ' 2 ) 

i 

< CAr 1 + CA- 2 . (4.3) 



Given t > 0, we choose A = i 1 / 3 and find P(n~ 3 / 2 £\ Wf > i) = 0(t~ 2 / 3 ). 

Hence, n -3 / 2 ^ VF 3 = O p (l), and it follows by (|4.ip and Markov's inequal- 
ity that X 4 = O p (l). 

Finally, we observe that, for any m > 4, 

P(u>(G(n, a)) > m) < p(x 4 > Q ) , (4.4) 

which thus can be made arbitrarily small (uniformly in n) by choosing m 
large enough. Hence, u>(G(n,a)) = O p (l). 

To complete the proof of Theorem [TJ^ii) let us note that for any fixed m < 
n, the probability that there are at least m vertices with weights W{ > n 1 / 2 
is larger than c\ > for some absolute constant c\ > 0, and conditioned on 
this event, the probability that the m first of these vertices form a clique is 
larger than C2 for some absolute constants c±, C2 not depending on n. 

Finally, we remark that all cliques of size four can clearly be found in time 
0(n 4 ). The number of such cliques is whp at most log log n, say, so there 
exists an algorithm which whp finds the largest clique in a polynomial time 
(for example by crudely checking all sets of cliques of size 4) . □ 
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5. a > 2 (finite second moment) 
Proof of Theorem Wiii) I Choose v such that 1/2 > v > I /a. Then (|2.2p (or 



2.1j) directly) implies that whp W max < n u . Furthermore, in analogy to 
421 and dOt . 



E^ttf; WW < ra") < E(^min(^,n^) 3 ) 

i i 

= n / 3x 2 P(VF > x)dx = 0(nn v ) = o(n 3/2 ), (5.1) 
J o 

and thus 

P ( n_3/2 IZ H/ i 3 > * 

i 

< f- 1 E(n- 3 / 2 Y W il ^max < n") + F(W max > n v ) 

i 

= o(l). (5.2) 
Hence, n~ 3 / 2 ^ ■ W 3 0, and it follows from (|4.ip that 

P(w(G(n,a)) > 4) = P(X 4 > 1) < E(min(l,E(X 4 | {Wi}?))) -» 0. 

Consequently, whp w(G(n,a)) < 3. 
Moreover, we can similarly estimate 

EX 3 <E Y Kj\kX jk = K{b 3 n- 3 Y W i W j W k) < ^ 3 (E^ 2 ) 3 ; 

i<j<k i<j<k 

(5-3) 

note that E W 2 < oo by (|2.ip and the assumption a > 2. Hence, the number 
of K% in G(n,a) is O p (l). To obtain the limit distribution, it is convenient 
to truncate the distribution, as we have done it in the previous section. We 
let A be a fixed large constant, and let be the number of K% in G(n, a) 
such that all three vertices have weights at most A, and let X%* be the 
number of the remaining triangles. Arguing as in (15. 3j) . we see easily that 

EX 3 A < i& 3 (E(Ty 2 ; W < A)) 3 (5.4) 

EI 3 A * < b 3 (E(W 2 )) 2 E(W 2 ; W > A). (5.5) 

Moreover, if W%,Wj < A, then Ay = 0(1 /n), and thus py ~ Ay, and it is 
easily seen that (|5.4p can be sharpened to 

El/ - ma = ^ 3 (E(Ty 2 ; W < A)) 3 . (5.6) 

Furthermore, we may calculate fractional moments E(X 3 ) m by the same 



method, and it follows easily by a standard argument (see, for instance, [13j, 
Theorem 3.19] for G(n,p)) that E(X 3 4 ) m — > /i™ for every m > 1, and thus 
by the method of moments 0, Corollary 6.8] 

X3 4 ^+ PoM (5.7) 
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as n — > oo, for every fixed A, where — — > denotes the convergence in distri- 
bution. 

Finally, we note that the right hand side of (15. 5ft can be made arbitrarily 
small by choosing A large enough, and hence 

lim limsupPpfg 4 * / 0) = 0, (5.8) 

A^oo n ^oo 

and that fiA — * H = g(&E(VF 2 )) 3 as A — ► oo. It follows by a standard 
argument (see Billingsley |3j, Theorem 4.2]) that we can let A — ► oo in (|5.7p 
and obtain 

X 3 ^Po(/i). (5.9) 
In particular, P(X 3 = 0) — > e~^, which yields the following result: 

¥(u(G(n,a)) = 2) -> e"" = e -|(^(w 2 )) 3 ; 

P(w(G(n, a)) = 3) - 1 - e"" = 1 - e ~KW 2 )) 3 . (5 - 10) 

Finally, note that G(n,a) whp contains cliques X2, i-e. edges, so clearly 
w(G(n,a))>2. □ 

6. Final remarks 

In this section we make some comments on other models of power-law 
random graphs as well as some remarks on possible variants of our results. 
We omit detailed proofs. 

Let us remark first that, for convenience and to facilitate comparisons 
with other papers, in the definition of G(n, a) we used two scale parameters 
a and b above, besides the important exponent a. By rescaling Wi 1— ► tWi 
for some fixed t > 0, we obtain the same G(n, a) for the parameters at a and 
bt~ 2 ; hence only the combination ab a l 2 matters, and we could fix either a 
or b as 1 without loss of generality. 

6.1. Algorithms based on degrees. As for the algorithmic result The- 
orem [2l it remains true if we search for large cliques examing the vertices 
one by one in order not by their weights but by their degrees and modify 
the definition of /C gr , /C qt , and /Cf t accordingly. (This holds both if we take 
the degrees in the multigraph G(n,a), or if we consider the corresponding 
simple graph.) The reason is that, for the vertices of large weight that we 
are interested in, the degrees are whp all almost proportinal to the weights, 
and thus the two orders do not differ very much. This enables us to find 
an almost maximal clique in polynomial time, even without knowing the 
weights. 

6.2. More general weight distributions. Observe that Theorems [Hand [2] 
remain true (and can be shown by basically the same argument), provided 
only that the power law holds asymptotically for large weights, i.e., (|2.ip 
may be relaxed to 

¥(W > x) ~ ax~ a as x -> 00. (6.1) 
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6.3. Deterministic weights. Instead of choosing weights independently 
according to the distribution W we may as well take a suitable deterministic 
sequence W{ of weights (as in Chung and Lu [7j ) , for example 

Wi = a x A* i = l,...,n. (6.2) 

All our results remain true also in this setting; in fact the proofs are slightly 
simpler for this model. A particularly interesting special case for this model 
(see Bollobas, Janson, and Riordan 0, Section 16.2] and Riordan [13]) is 
when a = 2, where (I2.3p and (j6.2H combine to yield 

ab 

Xij = 



/ij 

6.4. Poisson number of vertices. We may also let the number of vertices 
be random with a Poisson Po(n) distribution (as in, e.g., Norros fl3]). Then 
the set of weights can be regarded as a Poisson process on [0, oo) with 
intensity measure nd/i, where [i is the distribution of the random variable 
W in (|2.ip . Note that now n can be any positive real number. 

6.5. Different normalization. A slightly different power-law random graph 
model emerges when instead of (12. 3p we define the intensities Ajj by 

WiWj 

2^k=l W k 

(see for instance Chung and Lu |7[ and Norros and Reittu [15J]). Let us call 
this model G(n, a). In the case a > 1, when the mean E W < oo, the results 
for G(n, a) and G(n, a) are not much different. In fact, by the law of large 
numbers, ^^Wk/n — > Ely a.s., so we may for any e > couple G(n,a) 
constructed by this model with G(n, a)^ constructed as above, using (|2.3p 
with b = 1/(E W =F e), such that whp G(n, a)~ C G(n, a) C G(n, a) + , and 
it follows that we have the same asymptotic results as in our theorems if we 
let b = 1/EW. 

On the other hand, for a = 1, Yli ^fc = ( a + °p(l)) n l°g n ) an d fo r 
< a < 1, Y,i w k/n 1/a Y, where Y is a stable distribution with 
exponent a (e.g. see Feller 

EH Section XVII.5]). It follows, arguing as m 
Section El that for a = 2, the largest clique in G(n, a) has 



(l + Op(l))v / 2anlog" 1 n 

vertices, while for < a < 1 the size of the largest clique is always close to 
/n; more precisely, 



u(G(n, a)) 



Z = a2 Q / 2 y- Q / 2 , 



y/nlog a l 2 n 

where Z is an absolutely continuous random variable whose distribution has 
the entire positive real axis as support. (The square Z 2 has, apart from a 
scale factor, a Mittag-Lefner distribution with parameter a, see Bingham, 
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Goldie, Teugels 0, Section 8.0.5].) Thus, for a < 1, u(G(n, a)) is not 
sharply concentrated around its median; this is caused by the fact that the 
normalizing factor Wj is determined by its first terms which, clearly, are 
not sharply concentrated around their medians as well. Interestingly enough, 
since in the proof of Theorem [2] we dealt mostly with the probability space 
where we conditioned on W{, the analogue of Theorem [2] holds for this model 
as well. Thus, for instance, despite of the fact that neither the largest clique 
nor the full top clique are sharply concentrated in this model, one can show 
the sharp concentration result for the ratio of these two variables. 

6.6. The model min(Ay, 1). For small Ay, (|2,4p implies py ~ Ay. In most 
works on inhomogeneous random graphs, it does not matter whether we use 
(|2.4p or, for example, py = mm(Ay, 1) or py = Ay/ (1 + Ay) (as in Britton, 
Deijfen, and Martin-L6f ((J), see Bollobas, Janson, Riordan Q. For the 
cliques studied here, however, what matters is mainly the probabilities 
that are close to 1, and the precise size of 1 — p^ for them is important; 
thus it is important that we use (|2.4p (cf. Bianconi and Marsili [l|; 0] where 
a cutoff is introduced). For instance, a common version (see e.g. 0]) of 
G(n, a) replaces (|2.4p by 

Pij = min(Ay, 1). (6.4) 

This makes very little difference when Ay is small, which is the case for most 
i and j, and for many asymptotical properties the two versions are equivalent 
(see again [5]). In the case Yli Wf = o p (n 3//2 ), which in our case with Wi 
governed by (|2.ip holds for a > 2 as a consequence of __ fl5. 2j) , a strong general 
form of asymptotic equivalence is proved in Janson [12|; in the case a = 2, 
when Wf = o p (n 3 / 2 ) by (|4.2p . a somewhat weaker form of equivalence 
(known as contiguity) holds provided also, say, maxy- Ay < 0.9, see again 
[12|. In our case we do not need these general equivalence results; the proofs 
above for the cases a > 2 hold for this model too, so Theorem 
without changes. 

If a < 2, however, the results are different. In fact, (16. 4p implies that 
all vertices with W{ > b~ l l 2 n 1 / 2 are joined to each other, and thus form a 
clique; conversely, if we now define V~ = {i : W{ < (b + e) -1 / 2 ?! 1 / 2 }, then 
Pij = Xij < b/(b + e) for i,j G V~, and thus u(G(n, a)[V~]) = O(logn) whp 
by Lemma |U Consequently, arguing as in Section [3j 

u(G(n, a)) = (1 + o p (l))n¥{W > b^n 1 / 2 ) = (1 + o p (l))a6 a/2 n 1 - Q / 2 , 

so the logarithmic factor in Theorem [ ^i)| disappears. 

6.7. The model Ay/(1 + Ay). Another version of G(n,a) replaces (12.4p by 



(ii; nny hold 



Pij ~ l + Ay 



(6.5) 
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This version has the interesting feature that conditioned on the vertex de- 
grees, the distribution is uniform over all graphs with that degree sequence, 
see Britton, Deijfen, and Martin-L6f jf|. 

In this version, for large A^-, 1 — pij = 1/(1 + Xij) is considerably larger 
than for (|2.4ft (or (16. 4|) ). and as a consequence, the clique number is smaller. 
For a > 2, stochastic domination (or a repetition of the proofs above) shows 
that Theorem O p^rT pii) hold without changes. 



For a < 2, there is a significant difference. Arguing as in Section O we 
find that, for some constants c and C depending on a, b and a, whp 

cn (2- Q )/(2+a) < w ( G ( nia )) < Cn^^ 2+a \\ogn) a ^ 2+a \ 
Although this only determines the clique number up to a logarithmic factor, 

2- a 
2+a • 

2 ^ in Theorem [TJ 



note that the exponent of n is fr^, which is strictly less than the exponent 



6.8. Preferential attachment. Finally, let us observe that not all power- 
law random graph models contain large cliques. Indeed, one of the most 
popular types of models of such graphs are preferential attachment graphs 
in which the graph grows by acquiring new vertices, where each new vertex v 
is joined to some number k v of 'old' vertices according to some random rule 
(which usually depends on the structure of the graph we have constructed 
so far), see, for instance, Durrett @]. Clearly, such a graph on n vertices 
cannot have cliques larger than X n = m&x v < n k v + 1, and since for most 
of the models X n is bounded from above by an absolute constant or grows 
very slowly with n, typically the size of the largest clique in preferential 
attachment random graphs is small. 
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